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Abstract. - We drive a d-dimensional Heisenberg magnet using an anisotropic current. 
The continuum Langevin equation is analysed using a dynamical renormalization group and 
numerical simulations. We discover a rich steady-state phase diagram, including a critical point 
in a new nonequilibrium universality class, and a spatiotemporally chaotic phase. The latter 
may be 'controlled' in a robust manner to target spatially periodic steady states with helical 
order. 



How does an imposed steady current of heat or particles alter the dynamics of the isotropic 
magnet? To answer this question, we extend the equations of motion for the classical 0(3) 
Heisenberg model to include the effects of a uniform current in one spatial direction while 
retaining isotropy in the order parameter space. The resulting model is a natural generalization 
of the driven diffusive models of to the case of a 3-component axmZ- vector order parameter 
and, as such, is an important step in the exploration of dynamic universality classes Q far 
from equilibrium [Q. The form of the local molecular field in which spins precess in this 
driven state is strikingly different from that at equilibrium [P, and is responsible for all the 
remarkable phenomena we predict, including a novel nonequilibrium critical point and, in a 
certain parameter range, a type of turbulence. 

Here are our results in brief: (i) Despite 0(3) invariance in the order-parameter space, the 
dynamics does not conserve magnetization; (ii) As a temperature- like parameter is lowered, the 
paramagnetic phase of the model approaches a nonequilibrium critical point in a new dynamic 
universality class; (iii) Below this critical point, in mean-field theory without stochastic forcing, 
paramagnetism, ferromagnetism and helical order are all linearly unstable; (iv) Numerical 
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studies in space dimension d — I show spatiotemporal chaos in this last regime. This chaos, 
when 'controlled', is replaced by spatially periodic steady helical states. These predictions 
should be testable in experiments on isotropic magnets carrying a steady particle or heat 
current, as well as in simulations of a magnetized lattice-gas which we discuss at the end of 
this Letter. 

To construct our equations of motion, recall that at thermal equilibrium at temperature 
T the probability of spin configurations {Si} of a general nearest-neighbor Heisenberg chain 
with sites i is oc exp(— iJ/T), with an energy function 

H — — "^^JiSi.SiJ^i, (1) 

i 

where Ji is the exchange coupling between i and i + 1 . A spin Si s-t i precesses as Si = SiX hi 
where 

-> dH -> 

hi = — = JiSi+i + Ji-iSi-i (2) 

dS, 

is the local molecular field. Replacing Ji — > J{x) and Si — > S{x) in the continuum limit, yields 

§ 5^(2;) = J{x)SxdlS+{dJ/dx)SxdxS+.... For the physically reasonable case where J varies 

periodically about a mean value Jq, this reduces for long wavelengths to S{x) = JqS x d^S, 
which is invariant under x — > —x, even if the H is not. The dynamics conserves Si since 
it commutes with H . 

Now drive some background degrees of freedom in, say, the x direction, retaining isotropy 
in spin space. These could be some mobile species - particles, vacancies, heat - or some 
nonconserved internal variables. Possible microscopic realizations are discussed towards the 
end of the paper. For now, note that the dynamics in this nonequilibrium state does not 
follow from an energy function, and must be constructed anew. If we average over these 
background variables, their effect should be simply to modify the equations for the Si by 
allowing terms forbidden at thermal equilibrium. While many such terms are permitted, only 
two are relevant: (i) asymmetric exchange, i.e., hi = J+Si+i + J-Si-i yielding a precession 
rate gS x d^S + XS x dxS, with A oc J+ — J- proportional to the driving rate; and (ii) 
nonconserving damping and noise. Both (i) and (ii) were ruled out |^] at thermal equilibrium 
only because the dynamics had to be generated by (|l]) and (|^). Note that the A term, while 
rotation-invariant in spin space, is not the divergence of a current Q. The nonlinearity 
XS X dxS will thus generate nonconserving noise and damping terms even if these are 

not put in at the outset. 

For a general dimension d = dj_ + I, with anisotropic driving along one direction (||) only, 
the above arguments yield, to leading orders in a gradient expansion, the generalized Langevin 
equation 



dS 



d\\ + ri_yi ]S -vS - -{S ■ S)S - XS X diiS 



+g^^SxdlS + g^SxWlS + rf, (3) 

where we have allowed for spatial anisotropy in the coefficient of the usual spin precession term. 
The Gaussian, zero-mean nonconserving noise ff satisfies (77q,(x, t)77^(x', t')) = 2B SajS S'^'ix — 
x')6it - t') . 

In the equilibrium, isotropic limit, X = u = v = 0, = r± = r, the noise strength vanishes 
at zero wavenumber, and has a critical point where the renormalized r — > 0. In the driven 
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state, since the dynamics and noise are nonconserving, the critical point is w = 0, which in 
general takes place on a curve in the temperature/driving- force plane. As the drive is taken 
to zero there should be a crossover from nonequilibrium to equilibrium critical behavior. Our 
primary interest is in the behavior at a given nonzero driving rate, for which it suffices to 
vary the temperature-like parameter v in (^), keeping the rest fixed (with r||, r±,u > 0). For 
V > (the paramagnetic phase) all correlations clearly decay on finite lengthscales ~ ^/V^ 
and time-scales ^ 1/v, and nonlinearities are irrelevant. Let us focus first on the nature of 
correlations on the critical surface v = 0. Here, we expect an anisotropic scaling form for 
the correlation function C{x,t) = (S'(x x' , i -I- 1') • 5(x',t')): C(x,t) = x^'^F(i/a;|, a;_L/a;p), 
where i*" is a scaling function. In the linear approximation to (^) the roughening, growth, and 
anisotropy exponents are respectively x = 1 ^ d/2, z = 2, andC = 1, and F is analytic in its 
arguments. 

We now include the effect of the nonlinear terms in (^ via a standard implementation 
of the dynamical renormalization group (DRG) pO[ | based on a perturbation expansion in 
A and u. Rescaling x\\ — hxp x^ — &^x'^, t = h^t' and S = b^S', where b > 1 is an 
arbitrary parameter, the coefficients in Eq. (^) transform as r'^^ — fe^'^rn, = b^^'^'^r±, 
B' = b^-2x-Cid-i)-iB^ u' = b'^-'^u, a' = g = 6x+^-2g|| and 5^ = b^+^-'^^g^. A 

and u are thus relevant for dimension d < 4, and g\\ and g± are irrelevant for d near 4. In units 
where the ultraviolet cutoff is 1, A and u enter the perturbation theory in the dimensionless 
combinations r = {1/2it'^)\^ B / ^Jr^~r\ and k = {\/2tt'^)uB / ^ r^^r^, At the critical point 

z; = 0, setting the irrelevant g\\ and g±^ to zero, for d = 4 — e, we find to 0(e) the 

differential recursion relations 



I 

In ~ ' 4 



dl ^ ^ 48 

^ = S[z-2Y-C(d-1)-1 + — rl, 
dl ^ ^ ' 32 ^' 

Ot 35 

dn ,^ 11 , TT , 27 , n 

— = KiCe ttCk r) H -kCt . (4) 

Since we are working aXv — Q, we seek a fixed point that is stable with respect to perturbations 
in the remaining directions in the parameter space. For e = 4 — (i>Owe find the nontrivial 
stable fixed point r* = 64e/(357r), k* = 36 [1 + \/l409]e/3857r. The critical exponents for 
d < 4, to lowest order in e, are z = 2 — 16e/35, ^ = 1 — 2e/15 (anisotropic scahng) and [since 
u plays no role at 0(e)] X = 1 ~ d/2. These exponents clearly place this critical point in a 
new universality class. A more detailed analysis, including the approach to the critical point, 
will appear elsewhere 

We now investigate the low-temperature v < Q phase, in the absence 0/ noise. It is 
convenient to work with dimensionless variables, obtained by rescaling xj^, x\\, t and S in 
Eq. (^ ; this leaves A as the only parameter in the equation of motion. There are two 
static, spatially homogeneous steady states — a 'paramagnetic steady state' {Sa) = 0, and a 
'ferromagnetic steady state' (5*1) — {S2) = and {S3) — 1. It is straightforward to see from 
(^ that both these stationary solutions are linearly unstable We next look for static, 

spatially inhomogeneous steady states, a natural candidate being the helical state. Defining 
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Fig. 1 - Log- log plot of y = uj^\M:j(uj)\^ versus uj showing the l/uj'^ dependence of the power 
spectrum over approximately 1.5 decades. 



si and = tan-i(52/5i), 



for g\\ ^ g± ~ becomes 



v^^ + -iyp) ■ (V0) + -{SA\p - pd\\S^) , 



dt 
dt 



- V^S3 + S3-{p'' + Sl)S3 + Xp^d\\4>. 



(5) 



a, (f) — px|| and S3 — b (a, b and p are arbitrary constants) is a steady state 
-a2(l+A2)±^(a2(A2 + 1) - 1)2 - , and 2p = -A6± v^A^b^ - 4(i?2 - 1) 



A regular helix p = 
solution if 2fo2 = 1 

, where R = -i/a^ + W is the magnitude of each spin. The only free parameter a is bounded 
by a < (3 + A2)~-'^/2 fj-Qm the requirement that b be real. Unfortunately even this steady state 
shows a linear instability, triggered by the growth of S'3 @] . 

Having failed to find any stable static steady states analytically, we solve (||) numerically for 
d — 1 without noise, for a range of generic initial conditions. To avoid numerical instabilities 
we adopt an operator splitting method ||ll| — we solve the dissipative part using the standard 
Euler method and the drive part j|] by rotating each spin by an azimuthal angle |h(a;,i)|At 
about its computed local magnetic field h. With our choice of Ax = 1 and At = 0.0001 on a 
system of size N = 200 with periodic boundary conditions, we find that we avoid numerical 
instabilities and finite size effects. 

We find that the time series of the magnetization and energy density E = N^^ J dx{'VS)'^ 
never settle to a constant value; the motion could therefore be either temporally (quasi)periodic 
or chaotic. Figure.^ shows that the power spectrum of M3 = 5*3 goes as l/a;2. The power 
spectrum of E also shows a similar behavior. This suggests that the dynamics is temporally 
chaotic iQ. 

Space-time plots of local quantities, such as the signed local pitch, sgn(p) =sgn{dx(f>) (fig-§), 
strongly suggest the presence of spatiotemporal chaos . These results are preliminary, and 
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Fig. 2 - Space-time plot of the signed local pitch, sgn(p) = sgn{dx<j)) (black (+1), white (— 1)), 
revealing spatiotemporal chaos. 



only for d = 1. We shall characterize this behavior in greater detail elsewhere including 
studies of the dependence of the number of positive Lyapounov exponents on system size and 
the behavior for d > 1. 

The helix solutions of Eq. for v < are an infinite family of unstable spatially periodic 
steady states (parametrised by a) of the type discussed in jl^] . Can chaos in our model be 
controlled so as to stabilize and target p^ ] these helical states? The control of spatiotemporal 
chaos in PDEs is not nearly as well-developed as that in finite dimensional dynamical 

systems [|l3|. Accordingly, it is significant that we are able to stabilize, target, and hence 
control spatiotemporal chaos in our model, as we now show. 

For instance, in order to stabilize a specific helical configuration (with fixed a, b and p), 
we could in principle wait till the dynamics (presumably ergodic) eventually leads to this 
configuration, after which we apply small perturbations to prevent S3 from deviating from 
the value b (recall that the instability of the helical state was led by S3). This prescription 
successfully stabilizes the prescribed helix. 

In order to target this prescribed helix, we add to (^ terms which would arise from a 
uniaxial spin anisotropy energy V3 = r3(S'| — 6^)^ or V3 = r3{S3 — 6)^. We find that a 
sufficiently large and positive rs forces S3 to take the value 6 exponentially fast starting from 
arbitrary initial configurations. The subsequent evolution, given by Eq. (|^) on setting S3 = 5, 
can be recast as purely relaxational dynamics, 

dp__ 6F_ 90 _ ^_5F_ 

dt ~ ' at " ^ ^ 

where the 'free-energy functional' F has the form of a chiral XY model, 



F ^ 1 
2 



{vpr+p^iv4r-{p'+b') + ^-{p'+bY 



+A6p'a,|0] . (7) 
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Using the chain-rule, it is easy to see that dF/dt = \—{5F/5Si)'^ — {5F/6S2)'^\ < 0, Hence 

is a Lyapunov functional for the dynamics. Completing the squares, we see that 
appears in F in the combination (l/2)p^(9||(/) + A6/2)^, which is minimized by the helix 
(j) = — (l/2)A6a;|| . Starting from any initial configuration, the system plummets towards this 
unique helical minimum of F. 

Let us now see whether our control is robust against noise. We modify (|^) by including the 
noise rf in Eq. (^), and ask for the statistics of small fluctuations with Fourier components pk.{t) 
and (f>-k{t) about the controlled helical state, where 2tt/L < k < A for a system of linear extent 
L. It is clear from that the means (pk) and (0k) decay exponentially to zero: the relaxation 
time for pk is finite at small k, whereas that for ^k goes as k~'^. To calculate the variances, 
note that the dynamics is governed in the mean by the Lyapounov functional , and that the 
noise is spatiotemporally white. It follows that the steady-state configuration probability 
P[p,(j)] oc e~'^^ where c is an effective inverse temperature. F ~ /[const(p)^ -I- const(V(^)^] 
for small fluctuations about the helical minimum, i.e., P is approximately gaussian, so that 
(l^kP) ~ and (|pkP) ~ const, for small k. Thus the variance (fP) = /j.(|pkp) is L- 
independent for L — s- oo in any dimension d, whereas {(j)^) = /j,(|0kP) diverges as L and InL 
respectively for d — 1 and 2, and is finite for d > 2. Thus occasional excursions from the 
controlled state as a result of the noise do not lead to an instability of the targeted state for 
d > 2; the behavior for d < 2 is no worse than for a thermal equilibrium XY model. 

Having arrived at the continuum equations Eq. (j^) based only on symmetry arguments 
and conservation laws, we now suggest ways in which the driving nonlinearity in (^) may 
be realised, (a) Consider the isotropic magnet on a lattice whose unit cell lacks x — x 
symmetry. Now subject the spins to a spatiotemporally random, isotropic, nonconserving 
noise source. The lack of invariance under x-inversion together with this nonequilibrium noise 
should give rise, via an interplay of nonequilibrium noise and the electrons participating in 
the exchange interaction, to the S x dxS term. We admit, however, that we do not have a 
microscopic derivation of this, (b) Consider the asymmetric exclusion process (ASEP) |^ on 
a 1-dimensional lattice, each site i of which can be either vacant {ni{t) = 0) or occupied by 
one particle (ni(t) = 1) at time t. Each particle has an attached Heisenberg spin and may 
hop to the nearest neighbor at the right (left), if vacant, with probability p {q). The spin 
Si(t) at an occupied site i is the spin of the occupying particle. The local field at site i is 
hi{t) = Ji-i^iSi-i{t) + Ji+i^iSi+i{t). The required asymmetric exchange is achieved by making 
the Js dynamical. Recall that the exchange coupling between, say, i — 1 and i is operative 
at time t only if ni-i(t)ni{t) = 1. Imagine that the local exchange coupling emanating from 
a given site depended on some internal degree of freedom with a short relaxation time, as 
follows: Ji^i±i{t) = Ji if ni{t — 1) = 1; if not, Ji,i±i{t) ~ J2- Assume for simplicity that the 
particles can hop only to the right. Then a configuration 111 at sites i ~ 1, i, audi -I- 1 at 
time t was either already present at time t — 1 or arose from Oil (by a right hop from i — 2). 
Thus, for 111 configurations, Ji-i^i will be a weighted average of Ji and J2, while Ji+i^i will 
be Ji. In the continuum limit, and averaging over the particle dynamics, we will get the 
driving term in (||), with \ (x p — q. Such averaging is justified if the hopping species has a 
faster dynamics than the spins (for example, a smaller dynamical exponent). Alternatively, 
allowing evaporation-deposition in the ASEP renders the "particles" fast without altering 
qualitatively the derivation above for the effective asymmetric exchange. Of course, particle 
non-conservation induces spin non-conservation trivially in this case. The work of [| on 
moving space curves suggests another promising approach to finding realizations of our model. 

In conclusion, we have studied the interplay of dissipation, precession, and spatially 
anisotropic driving on the dynamics of a classical Heisenberg magnet in d space dimensions. 
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We have found a nonequilibrium critical point which we have shown, in an expansion in 
e = 4 — d, to be in a new dynamical universality class. We have presented evidence of spa- 
tiotemporal chaos in the mean- field dynamics of the model, at least in d = 1 , and have shown 
how this chaos can be controlled to yield helical order. Our work reinforces the idea that 
spatiotemporal chaos is a generic feature of driven, dissipative, spatially extended systems 
with nonlinear reactive terms. Further properties of this remarkable model, including the 
Lyapounov spectrum of the chaotic state, the possibility of complex ordered states or spa- 
tiotemporal chaos for d > 1, possible experimental realizations, and the crossover between 
equilibrium and driven behavior, will be discussed elsewhere 

* * * 
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